Friction and dilatancy in immersed granular matter. 
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The friction of a sliding plate on a thin immersed granular layer obeys Amonton-Coulomb law. 
We bring to the fore a large set of experimental results which indicate that, over a few decades of 
values, the effective dynamical friction-coefficient depends neither on the viscosity of the interstitial 
fluid nor on the size of beads in the sheared layer, which bears out the analogy with the solid-solid 
friction in a wide range of experimental parameters. We accurately determine the granular-layer 
dilatancy, which dependance on the grain size and slider velocity can be qualitatively accounted 
by considering the rheological behaviour of the whole slurry. However, additional results, obtained 
after modification of the grain surface by a chemical treatment, demonstrate that the theoretical 
description of the flow properties of granular matter, even immersed, requires the detailed properties 
of the grain surface to be taken into account. 

PACS: 47.57.Gc: Granular flow; 83.50.Ax: Steady shear flows, viscometric flow; 83.80. Hj: 
Suspensions, dispersions, pastes, slurries, colloids; 81.40.Pq: Friction, lubrication, and wear. 
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Conducting studies on immersed granular flows re- 
mains of primary interest. A host of geophysical or in- 
dustrial issues deals with mixtures of grains and fluid as 
submarine avalanches [H , snow flows or clay suspensions 
Q. Also fundamental issues are at stake: one would like 
to extend the empirical friction law proposed for dense 
and dry granular flows Q to immersed ones [1] . Enlight- 
ening previous studies of sheared and immersed granular 
media are numerous, and different devices have been de- 
veloped to describe mixtures of grains and fluids (mean- 
ing air or liquids). Among them we choose to focus on 
the three following. 

Studying immersed granular matter flowing down an 
inclined plane, C. Cassar et al. measured the dynamical 
friction-coefficient, for different flow configurations 
Their results were analyzed using an approach inspired 
by recent results obtained for dry and dense granular 
fiows 0, 01 : They report the friction coefficient as a func- 
tion of the dimensionless parameter /, first introduced 
by Da Cruz et al. Q, defined to be the ratio of an apt 
microscopic time scale (inertial, viscous, . . . ) to the rel- 
evant macroscopic time scale 7""'^, where 7 denotes the 
shear rate. For immersed granular- matter in the viscous 
regime I = {'i Vf)/{'^ Pg)^ where 77/ denotes the vis- 
cosity of the intersticial fluid, Pg the pressure exerted on 
the sheared media and a the normalized permeability of 
the granular packing [Ij]. They propose a semi-empirical 
law for which describes the whole set of data the 
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report for both aerial and immersed granular flows 

Using a Couette geometry Bocquet et al. tuned the 
pressure within the granular material by applying an up- 
ward air flow between the rotating and the stationnary 
cylinder [^. They found out that mean-flow properties 
and fluctuations in particle motion are coupled. They 
introduced an hydrodynamic model which quantitatively 
describes their experiments: The key feature is that the 
shear force obtained from this model is found to be pro- 



portionnal to pressure and approximately independent of 
the shear velocity. This model does not include any fric- 
tional forces between grains, but contains a phenomeno- 
logical relationship between the viscosity and the dilation 
of the media. 

Using an experimental setup flrst designed to perform 
sensitive and fast force-measurements in the dry case [9,] , 
Geminard et al. brought to the fore a dynamic friction- 
coefficient /i in the case of an immersed granular layer 
sheared by means of a sliding plate [l^ . At low imposed 
normal-stress, the friction force is shown to be indepen- 
dent of the plate velocity, which holds true as long as the 
granular material is allowed to dilate [11]. The main dif- 
ference with the dry case lays in the fact that the slider 
usually exhibits a continuous sliding instead of the stick- 
slip motion and in the value of the friction coefficient 
which is roughly cut down by half 13, The depen- 
dance of the frictional coefficient on the fluid viscosity 
and of the associated dependance of the dilatancy on the 
slider velocity were not reported. 

Here we report a set of experimental measurements 
of the friction coefficient and of the dilatancy in a wide 
range of fluid viscosities and grain sizes at very low I 
li]. Such a study is relevant for several reasons: First, 
we point out that the quasi-static regime is unaccessible 
to the free-surface-flow geometry as size effects crop up 
in this limit [4]. In addition, there is a strong discrep- 
ancy between the limit of fi for vanishing / reported in 
[H in ~ 0,43 with / ~ 4.10^3) and those reported for 
the plane-shear geometry by Geminard et al. {ji ~ 0.23, 
[13, (il) and S. Siavoshi et al. {n ~ 0.54, [l^), both for 
/ ~ 2.10"*. How can be explained such discrepancies 
between those three results? We also raise the following 
questions: What does happen when the fluid viscosity 
or the bead size are changed? How far does the anal- 
ogy with the Amonton-Coulomb laws remain relevant? 
What does the effective friction-coefficient depend on? 
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We choose to stick to the canonical plane-shear geom- 
etry for which we know that there is a strong analogy 
between the friction of a slider on an immersed granular 
layer and the Amonton-Coulomb law In the cho- 

sen geometry, the layer is free to dilate which makes it 
possible to measure both the friction coefhcient and the 
dilation of the granular layer at imposed normal stress. 

Horizontal 
Sensor Spring Vertical 
Sensor 




Granular 
Layer 



FIG. 1: Sketch of the experimental setup. 

Experimental setup. - The experimental setup (Fig.[T]) 
is very similar to the one described in [l^, A thin 
plate, the shder, is pushed at the free surface of an im- 
mersed granular layer by means of a steel leaf-spring 
(fc = 129 ± 2 N.m~^) connected to a translation stage 
driven at constant velocity, Vs by a computer-controlled 
stepping motor {Vg ranging from 0.1 to 100 /im.s^^). The 
coupling between the spring and the plate is insured by 
a metal bead, which avoids applying a torque. The fric- 
tional force is monitored by measuring the receding of 
the spring from its rest position with an inductive sen- 
sor (ElectroCorp, EMD1053). The dilatancy is obtained 
from the vertical displacement of the slider: A second 
inductive sensor, at rest in the laboratory referential, 
monitors the distance to a metallic target which endows 
the slider, which consists in a thin (5 mm) PMMA plate 
(76 X 53, 53 X 51, or 53 x 24 mm^). The granular ma- 
terial consists in spherical glass beads (Matrasur Corp.) 
sieved in order to obtain the three following mean di- 
ameters d = (100 ± 11), (215 ± 20) and (451 ± 40) ^lul, 
with a relative standard deviation almost independant of 
the characteristic grain size. The intersticial fluid con- 
sists in distilled water, water and sugar mixtures (viscos- 
ity r] ranging from 1 to 76 mPa.s), or Rhodorsil silicon 
oil (Rhodorsil, viscosity rj ranging from 71 mPa.s to 500 
mPa.s). All viscosities were measured, in addition, using 
an Ubbelohde viscosimeter. The thickness of the granu- 
lar bed (6.0 mm) is always larger than ten bead diameters 
in order to make sure that the sheared zone is not lim- 
ited by the bottom of the container and, thus, that edge 
effects are not at stake [l3| • Finally, the contact of the 
slider with the granular layer is insured by gluing a layer 
of the largest beads (451 ^m) onto the lower surface. We 



checked, for a layer of 215 /im-in-diameter beads, that 
the friction coefficient is independent of the size of the 
beads in the glued layer as long as it remains larger than 
that of the beads in the granular bed (table HI Top). 

(/im) 100 215 451 

fid ifim) 0.33 ± 0.02 0.38 ± 0.02 0.37 ± 0.02 
S (mm^) 53 x 24 53 x 51 76 x 53 
Aid (fim) 0.38 ± 0.02 0.42 ± 0.02 0.38 ± 0.02 



TABLE I: Top: Measured friction coefficient /i as a function 
of the diameter of the beads that insure the contact at the 
bottom surface of the slider (The sample consists of 215 fim 
beads in water); Bottom: Friction coefficient measured with 
sliders having different surface area and aspect ratio (100 fim 
beads in a water-sugar mixture, vicosity rj = 4.3 mPa.s.) 
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FIG. 2: Dynamic frictional force vs. effective mass of 
the slider. The effective mass of the slider is obtained by 
reducing the weight by the buoyancy force. From the slope 
one can infer /i = 0.41 ±0.02 (Vs = 3.5 fj,m.s~^ , 451 /xm beads 
in silicon oil, viscosity i] — 500 mPa.s). Inset: Friction coeffi- 
cient vs slider velocity Va in the same experimental conditions 
for the mass pointed by the rectangle. The size of the symbols 
indicates the error bars. 

Friction coefBcient. - In our experimental conditions, 
the value of the spring constant k is chosen so as to ob- 
serve the continuous motion of the slider in the whole 
accessible range of the driving velocity Vg . After a tran- 
sient regime, the frictional force reaches a steady-state 
value which is observed to scale up with the slider mass 
TO (Fig. [21), provided that the buoyancy force is taken into 
account, and to be independant of the slider velocity Vs 
(inset. Fig. [2]). In addition, we checked that n does not 
signihcantly depend on the slider surface-area or aspect- 
ratio (Table [H Bottom), as already known for the dry 
case [1J|. We repeated the procedure for different bead- 
diameter (d from 100 fim to 450 /im) and fluid viscosity 
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(?7 from 1 mPa.s to 500 mPa.s). We found out that fi 
neither depends on rj nor on d in the whole experimental 
range. In order to encompass those two results and the 
independence on the slider velocity Vs, we report ^ as a 
function of the Reynolds number Re = pdVs/r], where p 
stands for the fluid density. We estimate fi = 0.38 ±0.03 
for 10~^ < / < 5.10~^, which nicely supplement the data 
reported for the free-surface- flow configuration in [ij that 
limited to / > 4.10"^. 
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FIG. 3: Friction coefficient vs. the Reynolds number. 

We report data for three different diameters d: ■: 100 ^m; 
A: 215 ^m; □: 451 /^m (77 ranging from 1 to 500 mPa.s ). 
Note that fj, is constant in a range of Re covering more than 4 
orders of magnitude. Inset: Frictional force in the steady state 
regime vs slider position. The dynamic friction-coefficient is 
defined to be the mean value of the frictional force (14 = 
8.8/im.s~'^, rj = 1 mPa.s and m — 10.1 g). 



Layer dilation. - As we shall see, in contrast to the fric- 
tional coefficient, the dilation of the granular layer is sen- 
sitive to the grain size. Experimentally, the moot point 
consists in obtaining a reproducible reference state. The 
chosen procedure is as follows: In order to obtain a well- 
defined state of compaction, we initially push the slider 
over a distance of approximately 10 bead-diameters in 
the steady regime at a given velocity, henceforth denoted 
Vref (usually 8.8/im.s~^, excepted when specified). We 
then stop the translation stage and move it backwards 
until the spring goes back to its rest position without 
loosing contact with the slider, which remains at rest 
(contact loss could make the slider surf over the granu- 
lar layer as we push it forth at large velocity, meaning 
above AO iim..s~^) . We then immediately push the slider 
forwards at various driving velocities, Vs, over a few mil- 
limeters and monitor the vertical position of the plate. 
We observe that the total variation of the vertical posi- 
tion of the plate Ah, or total dilation, increases with the 
bead diameter d and the velocity Vg. By constrast, Ah 
does not significantly depend on the interstitial-fluid vis- 



cosity r] (Fig. |4]). We checked that these latter measure- 
ments neither depend on the preparation of the granular 
layer (by varying the velocity of reference Vref ) , nor on 
the slider mass m (Fig.[31 inset). 
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FIG. 4: Total dilation of the layer Ah vs. slider 
velocity Vs. (symbol, bead diameter, fluid viscosity): {+, 
100 /im, 1 mPa.s); (■, 215 ^m, 1 mPa.s); (A, 451 Atm, 1 
mPa.s); (□, 451 ^m, 71 mPa.s); (o, 451 ^im, 500 mPa.s); The 
dashed lines correspond to the interpolation of the experi- 
mental data to Eq. [T] with j3 = 2.9 ± 0.3. Inset: similar 
results for different reference velocity and normal stress: (■, 
100 /im, Vref = 3.5//m.s~^); (□, 100 ^m, Vref ~ 14/im.s~^); 
(o, 100 ^m, Vref = 3.5/im.s~^ and overloaded by 10.0 g). 

Discussion. - The dependence of the total dilation on 
the velocity and on the bead diameter can be accounted 
by the two following ingredients: First, we can guess that, 
due to the steric interaction between the grains (solid 
contact between grains or hydrodynamical interaction), 
the local shear-stress Cs induces a local normal-stress 
CT„ ~ Oi{(j)) as- In a first approximation, the coefficient a, 
which describes a geometrical property, depends only on 
the volume fraction of the grains, (p, and not on the shear 
rate, 7. The assumption is correct, at least for a dense 
suspension in the limit of small 7 [16[. We point out that, 
from this local relation between an and as , we recover the 
apparent friction law, F = pmg where n = 1/a, inde- 
pendant on the shear rate 7 or on the surface area S of 
the slider provided that a does not significantly depend 
on (/) Second, we assume, as already proposed by 

Bocquet ei aZ Q , that the rheological behavior of the im- 
mersed granular-material can be accounted, in addition, 
by as = ri((j>) 7, where the effective viscosity rj diverges 
algebraically as a function of near a critical volume 
fraction 0c: rj = '7o/(l ~ 4'/4'c)^ ■ At this point, assuming 
that the total dilation Ah is mainly due to the dilation 
of a constant number N of layers underneath the slider 
and linearizing the velocity profile in this region where 
the dilation is the larger, we write 7 = Vs/{Nd) and get 
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the following scaling law: 

Ah oc Vs'/^d^^'''>/^ 



(1) 



The interpolation of the experimental data reported in 
Fig. m leads to (3 = 2.9 ± 0.3, which compares quanti- 
tatively with the values reported in [§1 and references 
therein. 

Thus, our experimental measurements are in agree- 
ment with the conclusions of references [H and B. H ow- 
ever, the discrepancy with the values issued in [id . Il2 | 
and remains unexplained and deserves to be dis- 
cussed. Few physical origins can be taken responsible 
for such a dispersion related in the literature. Among 
them, in the case of spherical glass-beads, the polydis- 
persity of the batches and the properties of the beads 
surface, especially its roughness. On the one hand, using 
unsieved samples (polydispersity about 25%), we mea- 
sured significantly higher values n ~ 0.47 ± 0.02 of the 
friction coefhcient, which provides us with a rough esti- 
mate of the polydispersity effects. For all the data re- 
ported above, we use sieved samples, the standard devi- 
ation in the diameter being about 10% as in references 
13, 12, 13] J so that we estimate that the polydispersity is 
not enough for explaining the observed discrepancy. On 
the other hand, only a few observations, dealing with the 
influence of the surface properties on the effective fric- 
tion coefficient, have been reported: It has been recently 
proven that the roughness can drastically alter the dy- 
namic angles of repose for dry materials 15| and even 
be a motor for shear-induced segregation in immersed 
granular- materials [17|] . The beads used in our study and 



ties. In order to test the influence of the bead surface, 
we altered 451 /itm diameter beads by immersing them 
for 30 min in a 1.0 mol.L"^ sodium-carbonate solution 
[lit . The beads were then thoroughly washed up using 
distilled water and fi immediatly assessed. We obtained 
jj, = 0.30 ± 0.01 and, thus, a significant decrease of about 
20% by changing only the surface properties, even in im- 
mersed granular matter. We are currently performing an 
extensive and careful study of the friction coefficient as 
a function of surface roughness of the beads. 

As a conclusion, our measurements of the friction coef- 
ficient nicely supplement the results reported by Cassar 
et al [l| and we interpret our measurements of the di- 
latancy in the framework of the hydrodynamical model 
proposed by Bocquet et al. Q . The friction coefficient ^ 
remains constant over a large range of Re values because 
the dilatancy of the layer is a free parameter that adapts 
as the slurry is sheared at different velocities. In this 
regime, the effective friction-coefficient depends neither 



on the fluid viscosity nor on the bead diameter. How- 
ever, we underscore that the properties of the grain sur- 
face play an important role in the rheological properties 
of immersed granular-matter and could be responsible for 
the dispersion of /i values encountered in the literature. 
Finally, we emphasize that measurements of /i, defined 
from the mean value of the friction force in the steady- 
state regime, does not provide any piece of information 
about the grain- or fluid-characteristics. As an extension 
of this work, we are currently focusing on the fluctuations 
of the frictional force in the steady-state regime (inset. 
Fig. [3]) , from which we hope to recover a signature of the 
components of the slurry. 
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